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Despite tlie availability of a spin Hamiltonian for the Gd3Ga5 0i2 garnet (GGG) for over twenty 
five years, there has so far been little theoretical insight regarding the many unusual low temperature 
properties of GGG. Here we investigate GGG in zero magnetic field using mean-field theory. We 
reproduce the spin liquid-like correlations and, most importantly, explain the positions of the sharp 
peaks seen in powder neutron diffraction experiments. We show that it is crucial to treat accurately 
the long-range nature of the magnetic dipolar interactions to allow for a determination of the small 
exchange energy scales involved in the selection of the experimental ordering wave vector. Our results 
show that the incommensurate order in GGG is classical in nature, intrinsic to the microscopic spin 
Hamiltonian and not caused by weak disorder. 



The diversity of empirical data collected over the past 
fifteen years has demonstrated that geometrically frus- 
trated triangular and tetrahedral arrangements of an- 
tiferromagnetically coupled spins are highly partial to- 
wards the realization of exotic correlated phases in mag- 
netic materials [H 0, Q • The reason for the rich and typi- 
cally material specific properties of frustrated magnets is 
understood. It stems from their sensitivity to perturba- 
tions beyond the frustrating nearest-neighbor antiferro- 
magnetic (AFM) exchange which, on its own, leads to a 
macroscopic number of exactly degenerate and compet- 
ing, hence fragile, classical ground states. In this paper 
we show, through a careful theoretical analysis of neu- 
tron scattering experiments, that the extensively studied 
Gd3Ga50i2 garnet (GGG) is precisely such a system, 
though evidence for this fact emerges from a perspective 
on the problem that has heretofore escaped scrutiny. 

GGG displays a gamut of complex and interesting 
low temperature magnetic phenomena. In zero magnetic 
field, the behavior of GGG is uniquely rich. The nonlin- 
ear magnetic susceptibility xs peaks at Tg ^ 180 mK 
indicating a spin glass transition [^. However jinuon spin 
relaxation 6, 7] and Mossbauer spectroscopy [8| find per- 
sistent spin dynamics down to T <g; Tg. Meanwhile, pow- 
der neutron scattering data d, indicate that GGG is 
on the verge of developing true incommensurate long- 
range magnetic order with a correlation length (^ w 100 
A) extending over 8 cubic unit cells below 140 mK. 

A Hamiltonian Ti, describing GGG, that we shall ex- 
plicitly define below, has long been available ll|. It as- 
sumes classical Gd'^+ spins, is parameterized as a sum 
of empirical exchange contributions up to third nearest- 
neighbors as well as a magnetic dipolar contribution, and 
ignores potentially important quantum fiuctuations or 
disorder inherent to GGG Q . Previous numerical studies 
based on Ti have been unable to provide a quantitative 
explanation for the bulk [1| and dynamical [a[3,|8| prop- 
erties of GGG or the incommensurate spin-spin correla- 



tions that develop below 200 mK [l3| . This could be 
interpreted as evidence that exotic mechanisms involv- 
ing either quantum fiuctuations or disorder effects are at 
play in GGG. For example, it has been suggested that 
disorder, such as Gd^+ on Ga^+ sites, is responsible for 
nucleating regions with the spin order observed in neu- 
tron scattering in zero magnetic field [1, . 

In this letter we address the adequacy of the Hamilto- 
nian Ti. of Ref. |Tl| to explain the low temperature prop- 
erties of GGG. It has been demonstrated 0, [lH that H 
with only a nearest-neighbor exchange interaction can 
describe liquid-like spin correlations [3, [l3|. We show 
below that Ti. can also describe well the incommensurate 
spin-spin correlations and that there is no obvious reason 
to extend H. Rather, it is the correct procedure to treat 
its many competing small energy scales, in particular the 
long range dipolar interactions, and the unusually precise 
specification of weak second and third nearest-neighbor 
exchange couplings, that are foremost needed, and that 
have been missed so far. 

In order to correctly parameterize Ti and the resulting 
spin correlations we rely on neutron scattering experi- 
ments [H,[i31j which are the most direct and unambiguous 
probe of magnetic correlations. Our basic hypothesis is 
that conventional criticality controls the development of 
the sharp peaks observed in neutron experiments [lo| , 
and that, as in conventional magnets, a mean-field theory 
(MFT) treatment of the pertinent Hamiltonian should 
be able to capture the Gaussian regime in the approach 
to criticality. Our finding of a unique ordering wave 
vector qord shows that GGG, unlike the AFM nearest- 
neighbor Heisenberg pyrochlore, is not pathological and 
is amenable to such a soft-mode determination and cal- 
culation of the neutron scattering intensity profile. We 
show that GGG is very similar to the now well under- 
stood dipolar spin ice materials 0, where a MFT 
that takes into account the long-range nature of the dipo- 
lar interaction [l3| (i) identifies the correct ordering wave 
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FIG. 1: MFT J(Q) of GGG with Ja = -0.003 K, .h = 0.010 
K [O] and r = lO"" for = 3, 4, 5, 10, 1000. The Rc = oo 
profile is obtained using the Ewald method. The profiles are 
uniformly shifted to emphasize their strong dependence on Rc 
for small values, as well as their saturation to a well defined 
limit at large Rc- 



vector IJ], and (ii) correctly predicts the symmetry of 
the scattering close to the critical temperature. 

We examine the spin correlations in GGG by cal- 
culating the powder neutron scattering intensity I{Q) 
and comparing it with experimental data The 
scattering q-vector-dependent intensity is given by 
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and j span the sites of a garnet lattice (c.f 
in Ref. 
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of the vector separation between them, and enumer- 
ate distances between the sites starting from the near- 
est neighbors. The size of the GGG conventional cu- 
bic unit cell is a = 12.349 A [i]. The isotropic vectors 
Si have length \J S{S -\- 1) to classically approximate the 
S — 1 /2 Gd^+ spins. S^ are the spin components per- 
pendicular to q, -F(|q|) is the Gd'^"'" magnetic form-factor 
|1EJ |. and (. . .) denotes a thermal average. The spin cou- 
plings in TL include dipolar interactions TYdip of strength 
D, cut-off at a distance Rc, and exchange interactions of 
empirical strength J„, n < 3. We show below that Rc 
must be set to the true long-range limit _Rc = oo in or- 
der to obtain physically meaningful results. Throughout 
the paper, we keep D = 0.0457 K l3] and the nearest- 
neighbor exchange Ji = 0.107 K [ll| fixed. However, be- 
cause of the large uncertainties on J2 and J3 estimated 
from high-temperature analysis (—0.015 < J2 < 0.009 
K, -0.035 < J3 < 0.100 K) fnl, and the unreliability of 
the Monte Carlo procedure [Uj to refine them because of 
poor statistics (see also discussion in Ref. we con- 

sider these as free parameters to be adjusted by fitting 
the experimental I{Q) Q. 



We use a Gaussian MFT to obtain the soft modes 
and calculate I{Q) of the model H. Using a self-adaptive 
q-grid, we can calculate /(q) (in the thermodynamic 
limit) at an arbitrary positive dimensionless temperature 
T = T/T<!^PT - 1, where T^^'^ « 2 K and depends on 
J2 and J3. Due to the log |r| growth of the powder MFT 
Bragg peak intensities, as opposed to the l/r growth of 
the q-dependent intensities, r must be set rather small in 
order to work with the experimental profiles in the criti- 
cal regime (cf. Fig.[4|). We consider a finite q-space grid 
with 32'^ points in the first Brillouin zone, and construct 
a three dimensional cubic interpolation scheme to sepa- 
rately evaluate the numerator and the denominator of the 
MFT /(q) (cf. Eq. A37 in Ref. M). For q « q^rd + G, 
where Bragg peaks occur (G is a FCC reciprocal lat- 
tice vector), we back up the interpolation by calculating 
/(q) exactly [lB|. This ensures q-grid size independent 
results 0. I{Q) is computed by spherically averaging 
/(q) numerically. 

We find that 7i with dipolar interactions and exchange 
interactions beyond nearest-neighbor possesses a unique 
qord- Our numerical study suggests that the soft mode 
selection in the full TC is highly sensitive to the values 
of J2 and J3. This finding is in general agreement with 
results from Monte Carlo simulations of GGG [9','l05 (see 
also Ref. Uj for a comment on the sensitivity of the GGG 
specific heat on J3). However, our crucial new observa- 
tion is that qord is also very sensitive to an ad-hoc choice 
of the cut-off distance Rc in Tidip- We demonstrate this 
by plotting in Fig.[T]/((5) for various Rc for J2 = —0.003 
K and J3 = 0.010 K (values from For Rc < 0(10^), 

not only is qord a non- monotonous function of Rc, as the 
positions of the peaks indicate, but so is the unit cell 
magnetic form factor, as is evident from the significant 
changes in the overall shape of the diffuse scattering pat- 
terns near criticality (r <C 1). Fig.[l]demonstrates that Ti, 
with different low Rc gives completely different spin-spin 
correlations in the critical regime. Only for rather large 
cut-offs does qord become insensitive to Rc, approaching 
the Rc = 00 limit only when Rc ^ 10'^. Figure [T] sug- 
gests that it is essential to set Rc to its true Rc = 00 
limit to avoid inducing spurious ordered phases by treat- 
ing the lattice sum in Ti dip by a short cut-off method, as 
was done previously 0, [lO, [ll|, [l^, 17 1 . We henceforth 
acknowledge this and implement the Ewald method [lB| 
to process i?dip with Rc = 00. 

The scattering profiles in Fig. [1] and the q positions of 
the sharp peaks, even in the i?c = 00 limit, are incom- 
patible with the experimental data Q (cf. Fig. [H upper 
panel). However, the uncertainty in previously estimated 
values of J2 and J3 ll| allows fine-tuning them in or- 
der to obtain a better match with the low temperature 
data . To proceed, we are guided by the general under- 
standing of disordered systems near the critical bound- 
ary between long-range ordered and spin glass phases 0] . 
Specifically, we make the reasonable assumption that. 




FIG. 2: The optimization neutron-scattering penalty func- 
tion, V, in the J2 — J3 plane behaves discontinuously and 
shows minima in the black domains. Only for the couplings 
demarcated by a dotted trapezoid does the theoretical scat- 
tering profile reconstitute both the location of the sharp peaks 
and the overall shape of the broad and diffuse experimental 
scattering . The circle corresponds to the values of J2 

and J3 suggested in Ref. [ll|. The upper-left and bottom-right 
white regions correspond to different qord = structures. 

given the sizeable and saturated ^ ~ 100 A correlated 
regions below 140 mK the weak disorder in GGG 
freezes-in the q-dependence of quasi-critical incommen- 
surate correlated regions as Tg is crossed upon cooling. 
We therefore use a two-step fitting procedure to deter- 
mine J2 and J3. First, we establish whether the model 
Ti, allows for a qord compatible with the experiment. Sec- 
ond, we examine whether the model with a given { J2, J3} 
set is able to reproduce the overall structure of correla- 
tions both in the paramagnetic and "frozen-in critical" 
regimes Q. 

We find that qord for Rc = 00 belongs to the (hhl) re- 
ciprocal plane, with its location within the plane highly 
dependent on J2 and J3. To determine the optimum J2 
and J3, we generate a sequence of data points {Q} = 
\cioTd{J2, J3) + G|, and then match these points to the 
experimental sequence of sharp peaks in powder GGG 
data 0. To do so, we introduce a penalty function, 
^(•^2, J3), that provides a measure of mismatch between 
the first three clearly discernible and strongest experi- 
mental peaks at Ql = 0.64 1"^ = 0.85 A"^ and 
Q3 = 1.07 A~^ of Ref. I9I (see 43mK experimental data 
in top panel of Fig. d]) and the closest theoretically 
determined peaks: V = 100 x m.a.Xi=i ^2,3 ^Qi, where 
AQi — mingg|g-j \Q — Q*\ with the numerical factor 100 
chosen to set the scale, V = 1, for a maximum mismatch 
ofO.OlA-i. A map of V in the J2 — J3 plane is shown in 
Fig. [21 A good match between theory and experiment is 
denoted in black. The discrepancy between experiment 
and theory increases as the gradation of gray varies from 
black to white. 
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FIG. 3: Scattering profile for regions of the J2 — J3 plane 
where the positions Qi of the sharp peaks have been opti- 
mized; T = 1.6 X 10~^. Upper left and bottom right panels 
represent correspondingly located black domains in Fig. [2] 
The other two panels represent small domains in the central 
and bottom-left part of Fig. [S] The overall shape of I{Q) 
varies significantly, but all the profiles display sharp peaks 
agreeing with the experiment (shown by vertical lines). 



Figure[2]shows that even small (at the « 1% level of Ji) 
changes in J2 and J3 show up as large changes in the or- 
dering wave vector and, correspondingly, in V too. There 
exist only relatively small domains (shown in black) that 
minimize V. It turns out, however, that for each of these 
domains, the overall shape of the MFT powder scatter- 
ing profile changes noticeably even though qord remains 
almost constant at the value of 27r/a [0.29, 0.29, 0]. Simi- 
larly to the standard experimental procedure to analyze 
powder neutron scattering data, once qord is determined, 
one needs to solve for the magnetic structure. The sen- 
sitivity of the unit cell magnetic form factor to {J2t J3} 
leads to a significant qualitative change in the diffuse scat- 
tering pattern near criticality, hence offering a second 
optimization channel for the determination of { J2, J3}. 

Figure [3] shows the dependence of the theoretical I{Q) 
on { J2, J3} within the black domains of Fig. [2] A com- 
parison of such theoretical profiles with the experimental 
one reveals that only in the upper left black domain are 
the theoretical and experimental profiles similar; a strong 
discrepancy in other domains allows us to exclude them 
from further consideration. Within the upper left do- 
main, the match is visually better for rightmost values 
of J2 with the restriction that its value cannot cross the 
domain boundary at — 0.004K. 

We support and further refine the above analysis by 
fitting the MFT scattering profiles to the experimental 
I{Q) in the paramagnetic regime. For each value of J2 
and J3 we consider a x^-type penalty function for the 
overall shape of I{Q), which is adjusted by parameters 
responsible for a uniform and linear in \Q\ background 
terms, as well as by the overall scale factor. With t as 
a free fitting parameter, we do the fits in the interval 
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To conclude, we have found that, as in spin ices (ll.[l3l. 
an ad- hoc Rc cut-off of the dipolar interactions at 



FIG. 4: Positions of the first three strongest peaks of the ex- 
perimental scattering profile (pointed to by arrows in the up- 
per panel) as well as the overall shape of scattering profile are 
used as criteria for the optimization procedure in the J2 — J3 
plane. The resulting mean-field theory profile at dimension- 
less temperature r — f .6 x fO~^ (lower panel) compares well 
with the experimental one. The positions of three experi- 
mental peaks are denoted by vertical lines. The experimental 
scattering profile has been obtained by digitizing Fig. 2 of 
Ref. H;. The theoretical profile is shown for J2 = — 0.005K, 
J3 = O.OIOK, denoted by a box in Fig. [1 



less than a few hundred nearest-neighbor distances leads 
to spurious long-range ordered phases. Only a proper 
treatment of an infinite Rc allows one to get a handle 
on the very small exchange interactions beyond near- 
est neighbors which dictate the incommensurate ordering 
wave vector qord • With a reasonably well parameterized 
Hamiltonian now in hand, further theoretical studies to 
explore the origins of the complex phenomena displayed 
by GGG in zero and nonzero magnetic field may now be 
possible. Perhaps most importantly, we have identified 
the likely reciprocal plane that contains Qordi with qord ~ 
27r/a [0.29, 0.29, 0]. With this prediction available, it may 
now be possible to perform single crystal neutron scat- 
tering measurements in a reflection geometry to beat the 
^^^Gd absorption problem in hig h quality single crystals 
with natural Gd abundance [l8|, and allow for a quan- 
titative investigation of the development of correlations 
in GGG at low temperatures. This may prove a fruitful 
endeavor to help shed light on the mysteries of GGG. 
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Q = 0.26 — 1.57A~^ which allows for the complete shape 
determination of the first experimental broad diffusive 
scattering peak at 1.05A~^ without a jQp background 
correction. Most importantly, an analysis of the (para- 
magnetic) data at 175 mK ^] allows us to confirm the 
outcome of the analysis of the critical data at 43 mK and 
reject all but the aforementioned upper left domain of 
Fig. [21 Within that domain the data impose conserva- 
tive error bars on J2: -0.012 K < J2 < -0.004 K, with 
limits on J3 set by that same domain, i.e. —0.003 K < 
J3 < 0.012 K, see Fig. H 

The MFT I{Q) with the optimized Jn values captures 
the features of the experimental I{Q) well (see bottom 
panel of Fig. 3]). The optimization procedure to match 
the three strongest experimental peaks leads to a theo- 
retical sequence of peaks that can be identified on the ex- 
perimental profile as rather well-distinguished peaks and 
cusps. The match (Fig.[4|) between MFT and experimen- 
tal I{Q) suggests post-factum that the qord dependence 
of the correlations can be described by MFT. 
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